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ABSTRACT 

Using ACS-HST images to yield continuum subtracted photometric maps in Ha of the Sbc galaxy 
M51 and the dwarf irregular galaxy NGC 4449, we produced extensive (over 2000 regions for M51, 
over 200 regions for NGC4449) catalogues of parameters of their Hii regions: their Ha luminosities, 
equivalent radii and coordinates with respect to the galaxy centers. From these data we derived, for 
each region, its mean luminosity weighted electron density, (rig), determined from the Ha luminosity 
and the radius, i?, of the region. Plotting these densities against the radii of the regions wc find 
excellent fits for {rie) varying as This relatively simple relation has not, as far as we know, 

been predicted from models of Hii region structure, and should be useful in constraining future models. 
Plotting the densities against the galactocentric radii, r, of the regions we find good exponential fits, 
with scale lengths of close to 10 kpc for both galaxies. These values are comparable to the scale lengths 
of the Hi column densities for both galaxies, although their optical structures, related to their stellar 
components are very different. This result indicates that to a first approximation the Hii regions can 
be considered in pressure equilibrium with their surroundings. We also plot the electron density of 
the Hii regions across the spiral arms of M51, showing an envelope which peaks along the ridge lines 
of the arms. 

Subject headings: ISM: general — HII regions — galaxies: structure 



1. INTRODUCTION 

Electron densities for Hii regions can be measured in 
two different ways, which respectively yield two differ- 
ent ranges of values. Density sensitive emission line 
ratios, such as the [SII] AA6731/6716A doublet ratio, 
or the [Oil] AA3729/3 726A doublet ratio give values in 
the range of hundreds fZarit skv et~aI1ll994l : lMcCan et al.l 
[198 5i ) . However measurements using Ha surface bright- 
ness to derive the emission measure, together with the 
estimated ra dius of an Hil regio n, give values in the range 
1-10 (see e.g. iRozas et al.lll996l'). This type of dichotomy 
was recognized long ago bv lOsterbrock fc Flathed (|1959l ^ 
who used the emission line ratio of the [Oil] doublet, 
AA3729/3726A, as their version of the first technique, 
and cm wavelength radio emission (from the literature) 
to estimate the emission measure, i.e. for the second 
technique. They reached the conclusion that the best 
way to reconcile the differences, of more than an order of 
magnitude between the two techniques, is to assume that 
the emission nebula is clumpy, with dense clouds embed- 
ded in a much more tenuous substrate. The line ratios are 
strongly weighted by density, and yield values reflecting 
the densities of the dense clumps, while the electron den- 
sity derived using the emission measure and the overall 
size of the Hii region is a geometrical average, whose bal- 
ance between the contribution of the clumps and the sub- 
strate depends on the fractional volume occupied by the 
former. In models for Hii regions based on recognition 
of this inhomogeneity, the fractional volume occupied by 
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the clumps is termed the filling factor, and a first or- 
der approximation can be made that the interclump sub- 
strate makes a negligible contribution to the line emis- 
sion, so that the interclump volume can be considered to 
have negligible electron density (see Ostcrbrock (198^) 
for a standard treatment). In the present article we are 
using the measured luminosities of Hii regions in Ha to 
compute an electron density using the second method, 
so our values are to be taken as the means produced by 
the second technique. We will use them to derive com- 
parative electron densities of selected Hii regions in M51 
and NGC 4449, to establish functional relations between 
them, the sizes of the Hii regions, and their positions 
within the discs of their respective galaxies, and hence 
to provide an overview of the global behavior of the elec- 
tron density within the discs of the two galaxies. 

2. THE OBSERVATIONAL DATA AND ITS 
TREATMENT 

M51 was observed in January 2005 with the ACS 
on HST using the broad band filters F435W (close to 
the standard B band) F555W (close to V) and F814W 
(close to I) and the narrow band filter F658N (Ha) 
within observational program 10452 (P.I. S. Beckwith, 
Hubble Heritage Team). The data used here were 
corrected, cali brated, and combined into a mosaic of 
the galaxy by iMutchler et~al] (|2005l ). The images of 
NGC 4449 were also taken with the ACS, through the 
same filters as those used for M51, in November 2005, 
within program GO 10585 (P.I. A. Aloisi). The im- 
ages were taken at two different positions along the ma- 
jor axis of the galaxy, with four d ithered pointings on 
each position (|Annibali et al.l [20081 ) to help cosmic ray 
removal. We c ombined them using MULTIDRIZZLE 
(jKoekemoer et al...2002 ). yielding overall exposure times 
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of 3600s, 2400s, 2000s and 360s respectively for the four 
named filters. The total resulting field is sized 345x200 
arcsec^; we produced mosaics containing complete im- 
ages of the galaxy in all bands. 

For both galaxies standard correction and calibration 
procedures were used (bias, flat-field, dark current, and 
distortion correction s) using the AC S pipeline calibra- 
tion suite CALACS (|Hack et alll2000h . Use of the Mul- 
tidrizzle process allowed us to correct for detector de- 
fects and distortions, though we retained the original 
pixel scales (0.05 arcsec per pixel). After this process 
we converted the images from e-/s to flux in erg cm~^ 
s~^ A "-^ using the standar d PHOTOFLAM conversion 
factor (jSirianni et al.l [20051 ). We aligned the continuum 
and Ha images and produced a continuum subtracted 
Ha image, after deriving the factor of proportionality 
between the continuum in t he on-band and off -band fil- 
ters, following the method of lBoker et al.l (| 19991 ). A more 
detailed description of our image treatment can be found 
in Gutierrez et al. (2010, in preparation). The method 
for measuring the Ha flux of an Hii region requires defln- 
ing its boundaries in the presence of any local diffuse Ha 
emission and any overlapping neighbour regions, using 
circular or rectangular deflning apertures. The bound- 
ary of an HH region is deflned as containing only pixels 
whose brightness is greater than 3cr, measured on the 
background. For NGC 4449 the diffuse background sub- 
traction presented special difficulties due to its strength, 
and to perform this accurately we used a selective fil- 
ter, as defined in Gutierrez et al. (2010, in preparation). 
Once a local background brightness was determined, it 
was multiplied by the number of pixels within the Hii 
region and subtracted from the flux directly measured 
within the region, to yield the measured flux from the 
region. To allow for the overlap of regions, if the bright- 
ness between the two maxima falls below 2/3 of the value 
of the fainter maximum we considered them to be two 
separate regions, with the boundary in the depression of 
the brightness distribution. Where this criterion was not 
met we classified the region as single. 

To allow for the contribution of the [Nil] doublet 
through the F658N filter, for M51 we used the mean 
of the observed ratios of these lines to Ha for 10 Hii 
regions in that galaxy by 'Bres olin et al.l (|2004f ). which 
show very limited variations from region to region, to 
calculate the factor by which the observed fiux should 
be reduced, using the filter transmission curve and the 
observed redshift ofM51. Fo r NGC 4449 we used the ob- 
served ratio from iKennicuttI ()1992l ). The factor for M51 
was 0.67, and the corresponding factor for NGC 4449 is 
0.86. We understand that radial decline in metallicity 
(especially in M51), will give rise to errors in using these 
blanket numbers for each object, but we can show that 
these errors are of second order (Gutierrez et al. 2010, 
in preparation). Also, we estimated the possible system- 
atic error due to the presence of the [OIII] emisison line 
at A5007A in the F555W filter and found this to be less 
than 1%. 

Finally, in deriving absolute fluxes we took the dis- 
tance to M51 as 8.39 Mpc, following iFeldmeier et al.l 
(|1997D and the dist ance to NGC 4449 as 3.82 Mpc, fol- 
lowing lAnnibair^Fal. (2008). We measured the luminos- 
ity of each region as outlined above, and also, from the 



continuum- and background-subtracted Ha image, de- 
termined the area subtended by the image, and derived 
a value for an equivalent radius, obtained dividing the 
area by tt and taking the square root. With these mea- 
surements we could determine a mean electron density 
for a each region, as described below. A catalogue of 
2657 regions of M51, and 273 regions of NGC4449 was 
produced, containing positions, measured fluxes and lu- 
minosities in Ha, and equivalent radii (see Gutierrez et 
al. 2010 for detafls). 

3. THE MEAN ELECTRON DENSITY OF THE HII 
REGIONS 

3.1. Electron density as a function of HII region radius 

We described in the introduction the inhomogeneity of 
Hii regions, but with the set of observations described 
here we conflne our attention to the mean el ectron den- 
sity ((rte)), and use the "Case B" formula (jOsterbrockl 
[19891 ) 
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to derive this, where L(Ha) is the Ha luminosity, a^ is 
the case B recombination coefficient, hvHa is the energy 
of an Ha photon, and R is the equivalent radius of the 
Hii region. If we measure R in pc and L(Ha) in erg 
s^^, and taking a canonical value for the temperature of 
10000 K we can rewrite equation [1] as 



(ne) = 1.5 X 10 
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In Fig. [T] we show the mean electron densities of all the 
regions observed in M51 against their equivalent radii. 
There is considerable scatter, which we will investigate 
further below. A linear fit to the points gives the result: 
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in qualitative terms, the larger regions have lower mean 
electron densities. The functional trend is quite well re- 
produced by the Hii regions in NGC 4449 as we can see 
in Fig. [21 With far fewer points, the best linear fit is 
given by 
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These are only two galaxies, but the trends can at this 
stage be usefully summarized by the simple expression 
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This relation is obeyed more closely for both galaxies 
if we correct the relationship between (rie) and R for 
the dependence of (rie) on the galactocentric radius, r, 
as described in section 3.2 below. Making this correction 
gives adjusted values of -0.52, and -0.47 for the exponents 
in M51 and NGC 4449, respectively. 



Distributions of mean electron density in HII regions. 
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Fig. 1. — Mean electron densities of all the Hll regions measured 
in M51 as a function of equivalent radius. In the lower panel the 
densities have been corrected by a factor which takes into account 
their galactocentric distances (plotted in Fig. [Sj. 
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Fig. 2. — Mean electron densities of the Hi regions measured in 
NGC 4449 as a function of equivalent radius. In the lower panel the 
densities have been corrected by a factor which takes into account 
their galactocentric distances (plotted in Fig. |^. 



3.2. Electron density as a function of galactocentric 
distance 

In the upper panel of Fig. [3] we have plotted the mean 
electron density, (ng), against galactocentric radius, r, 
for M51, including all the regions, but distinguishing by 
point styles among those found in different large scale 
features of the galaxy: arms, interarm zones, and the 
central kpc. The outstanding feature of this plot is the 
single straight line fit (which is an exponential, as the 
electron density scale is logarithmic) , which has the form: 

K) = K)„ e-^/\ (6) 

where {n^)^, with value 10±1 cm^'' is a central value of 
electron density and /i is a scale length, which takes the 
value lOibl.O kpc. We have limited the fit to galactocen- 
tric radii within r — 10.4 kpc, because beyond this radius 
there is an abrupt local rise and fall in (rig), attributable 
to the effect of the interacting neighbour galaxy NGC 
5195. We can make separate estimates for the constants 
(ue)^ and h in the arms and in the interarm zones. For 
the arms this gives us {n^)^ = 10 cm^'^ and h = 11 kpc, 
while for the interarm zones we find (rie)o ^ 8 cm^"^, and 
h = 11 kpc . 

Since we know that in the upper panel of Fig. [3] much 
of the apparent scatter is due to the (inverse) dependence 
of (rie) on the radii of the individual Hii regions, we next 
opted to apply a normalized correction factor for this ef- 
fect to all the data, and replot the figure, yielding the 
central panel of Fig. [3l in which the (exponential) linear 
fit is clearly much improved. To make the fit we have 
chosen to leave out the zones from 1.4 to 4.6 kpc from 
the center, and from 10.4 to 11.5 kpc from the center, 
and now find values of {ne)^ and h of 12 cm~^ and 9 kpc 
respectively. It is notable that the scale length found 
for the neutral atomic hydrogen component of M51 by 



iTilanus fc AllenI (|1991l l is 9.1 kpc, on which we will com- 
ment further below. We should point out that the zone 
between the center of the galaxy and a radius of 1.4 kpc 
does fit the exponential for the disc as a whole. 

To take a look at the global behavior we have taken 
averages of all the regions within annuli of width 200 pc 
and have plotted, in the lower panel of Fig. [3] these val- 
ues against galactocentric radius. This shows the overall 
trend in (ng), and we can pick out the general decline 
with radius, but note a plateau between 1.4 kpc and 4.6 
kpc. If we omit the points in this interval from those 
entering the exponential fit, we obtain values for (?^e)o 
and h, of 11 cm~'^ and 9 kpc respectively. 

To compute the effect of the temperature gradient 
on the determ ination of (rie) and of its gradient, we 
used data from lBresolin et all ()2004f ) giving the O abun- 
dance gradient, and the widely used CLOUDY model 
(iFerland et al.l[l998l ) to derive the effective temperature 
gradient, which we then incorporated into the expression 

(n,) = 2.3 X 10-i«y (7) 

and recalculated the electron densities. This introduced 
a change in the exponent of the relation between (rie) 
and radius from -0.56 to -0.55, and in the scale length 
from 9.2 to 9.9 kpc. These are the values represented in 
Fig.H 

Here we should point out the similarity of the lower 
panel of Fig. |3] with the radi al plot of the azimutha lly 
averaged Hi column density in ITilanus &: AllenI ()1991l ). 

General expressions for (rig) are then: 
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Fig. 3. — Mean electron density ((rie)) vs galactocentric radius, 
r, for M51, distinguishing by the point styles those found in the 
arms, in the interarm zones, and within the central kpc. The points 
in the centre panel were found by applying a normalized correction 
factor to those in the upper panel, in order to take into account the 
variation of (rie) with region radius, R. In the lower panel we show 
the result of taking the median value of (rte) for the full sample of 
regiones within annuli of width 200 pc. 
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Fig. 4. — Mean electron densities ((rie)) for the measured Hi 
regions NGC 4449. In the upper panel (rie) is plotted against 
galactocentric radius, in the centre panel (rie) has been modified 
by a correction factor to take into account the variation with the 
radii R of the individual regions (a l/R^^'^ dependence), and in the 
lower panel we plot the median values of all the regions sampled 
within annuli of 100 pc width. 



{45.8 e-*^/!" R-"-^^cm-^- 
r < 1.4 kpc & r > 4.6 kpc 
28.9 i?-0-55cm-3; 
r e [1.4,4.6] kpc. 

The results for NGC 4449 are summarized in Fig. 01 
In the upper panel of Fig. |3]we have plotted (rie) against 
galactocentric radius for all our measured regions, in the 
central panel we present the data normalized to take into 
account the l/i?"-^ relation explained above, while in the 



lower panel we have averaged the points in the central 
panel in annuli of width 200 pc. In the latter two panels 
the linear fits show the exponential radial decline in mean 
electron density, which is fitted by the whole galaxy out- 
side the central 120 pc. Using the best fit in the lower 
panel of Fig. |3]and normalizing, (rig) is given by: 

(ne) = 45.7 e-''/" ii'-o-^^ (9) 

This scale length here is 11 kpc which is muc h bigger than 
the op tical radius of the galaxy. However, iHunter et all 
(j 19981 ) observed an Hi disc with a central dense compo- 
nent some 4.5 kpc in radius embedded in a much larger 
elliptical component with a major axis limiting radius 
of 18 kpc. Taken our cue from M51, which has an Hi 
scale length of 9 kpc a nd a limiting Hi radius of 15 kpc 
(jMeijerink et al.ll2005l ). we can claim here that the Hii 
region scale length and the Hi scale length of NGC 4449 
are comparable, consistent with the scenario where the 
electron densities in the Hii regions follow the column 
density of the Hi gas, as in M51. 

3.3. The mean electron density as a function of 
position in the spiral arms of M51 

Finally in this section we show in Fig. [5] the electron 
density in the spiral arms of M51 as a function of the po- 
sition relative to the central ridge line of the arms. The 
ridge lines of the arms, used to define this relative posi- 
tion, were measured using the continuum image through 
the F814W filter. The abscissa in these plots is the dis- 
tance from the ridge line to the regions. We can see 
that the mean electron density reaches peak values along 
the ridge line, which is not surprising, but taking the 
complete populations of Hii regions we can see that the 
systematic increase of density towards the ridge line is a 
property of the upper envelopes of the plots for the two 
arms. There are many regions close to the ridge line with 
relatively low densities, though no high density regions 
towards the edges of the arms. 

4. CONCLUSIONS 

We have used the Ha luminosities and measured radii 
of the populations of Hii regions in M51 and NGC 4449, 
from HST images, to derive their mean electron densi- 
ties, aware that their inhomogeneous structures imply 
that the mean electron density does not describe either 
the in situ densities of their strongly emitting clumps, or 
the densities of their more tenuous interclump gas, but 
is a luminosity weighted mean over the Hii region vol- 
ume. We have found two functional dependences of this 
mean electron density parameter, valid for both galax- 
ies: it varies inversely proportionally to the square root 
of the radius of an Hii region, and it falls exponentially 
with galactocentric radius. The former relation must de- 
pend on the effects of the energy outflow from the mas- 
sive ionizing stars on their environment. As far as we 
are aware it does not correspond to the predictions of 
any p reviously published specific model (jScoville et al.l 
(|2001l ) show a graph for M51 which could be used to 
infer a similar result, but did not derive an explicit rela- 
tion for these variables), and sets interesting constraints 
on the physics of OB cluster formation and the subse- 
quent interaction with the cluster environment. The lat- 
ter relation, with the scale length for the electron density 



Distributions of mean electron density in HII regions. 
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Fig. 5. — The mean electron densities of the HII regions in M51 plotted against distance from the central ridge line of the arms in M51 
(squares). Also shown is the variation of the average value in 200 pc bins of these densities as a function of this distance (thick lines) for 
both arms. 



comparable with that for the Hi column density in galax- 
ies of very different sizes and types (M51 is a Sbc while 
NGC 4449 is an IBm), demonstrates that, globally, an 
Hii region is in quasi-pressure equilibrium with its sur- 
rounding gas. The quantitative differences between the 
density of arm and interarm regions are not great, but 
also reflect this quasi-equilibrium, since the values are 
systematically somewhat lower for the interarm regions. 



Finally the distribution of mean electron densities across 
a spiral arm has an upper envelope peaking along the 
center, the ridge line of the arm, but values well below 
the envelope are distributed fairly uniformly across the 
arm. These relations should be of value for improving 
physical descriptions of Hii regions as zones of interac- 
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